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Abstract 



In this article it is shown that in an equiUbrium classical canonical ensemble of 
molecules with two-body interaction and external field full Gibbs distribution can 
be uniquely expressed in terms of a reduced two-particle distribution function. This 
means that while a number of particles N and a volume V are fixed the reduced 
two-particle distribution function contains as much information about the equilib- 
rium system as the whole canonical distribution. The latter is represented as an 
absolutely convergent power series relative to the reduced two-particle distribution 
function. As an example a linear term of this expansion is calculated. It is also 
shown that Gibbs distribution function can de expressed in terms of reduced distri- 
bution function of the first order and pair correlation function. That is the later two 
functions contain the whole information about system under consideration. 

1 Introduction 

In classical statistical mechanics an equilibrium system of molecules in a volume V 
is described by canonical distribution function Fisf{q,p), where {q,p) is a set of phase 
variables: coordinates qi and momenta Pi of molecules. If interaction of molecules is 
additive, reduced distribution functions are introduced [lll2]- They are used for evaluation 
of thermodynamic characteristics of this molecular system. It's usually accepted that 
reduced distribution functions contain information about a molecular system less than 
the initial canonical distribution function. It's also supposed that the lower an order of a 
reduced distribution function is, the less information it contains. But there does not exist 
a proof of this statement in scientific literature. 

On the other hand, it is known that for an equilibrium canonical ensemble of non- 
interacting particles a canonical distribution function F/v(g,p) is decomposed into a prod- 
uct of reduced one-particle distribution functions Fi{q,p) [U [2]. This means that all 
information about such system is contained in the reduced one-particle distribution func- 
tion. 

In the article p] it was proved that for a system with pair interaction and without 
external field a reduced two-particle distribution function defined in terms of nonormal- 
ized Gibbs distribution (without configurational integral) contains the whole information 
about such system. In ^ that result was obtained for normalized Gibbs distribution and 
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respective reduced distribution function. This results were obtained for a system which 
is not subjected to influence of external field. 

In this paper it's proved that for such system placed into external field one-to-one 
correspondence between full canonical distribution function and a reduced two-particle 
distribution function exists. Moreover it's shown that one-to-one correspondence between 
a full canonical Gibbs distribution and a set of reduced distribution function of the first 
order together with pair correlation function exists too. This means that the reduced 
two-particle distribution function as well as mentioned set contain the total information 
about system under consideration. 

We consider an equilibrium system of particles contained in the volume V under 
the temperature T. Potential energy of system is supposed to have the form 

N 

f/7v(gi,...,g7v) = X^Mi(gi) + MQj^Qk), (1) 

i=l ^<j<k<N 

where ui{q) is an external field and U2{q, q') is a two-body interaction energy of particles. 
Probability distribution function of equilibrium system is the canonical Gibbs distribu- 
tion which is decomposed into a product of a momentum distribution function and a 
configurational one [H [2] . The former is expressed as a product of one-particle Maxwell 
distributions, the latter has the form 

DN{qi, ...,qN) = Qn^ exp{-/?t/iv(gi, . . . , giv)}, (2) 
where j3 = 1/kT, k is the Boltzmann constant and Qat is the configuration integral 

Qn = Jexp{-pUN{qi, ■ ■ . ,g7v)}rfgi ■ --dqN- (3) 

Here and below integrating with respect to every configurational variable is carried out 
over the volume V. For a system having interaction of form ([1]) reduced distribution 
functions are introduced by expressions [2] 

N\ f 

Fi{qi, ...,qi)= _ J DN{qi, • • • , qN)dqi+i ■ ■ ■ dq^, I = 1,2,.... (4) 

These functions are used instead of full canonical distribution ([2]) to calculate various 
characteristics of the molecular system. Let us investigate properties of the reduced two- 
particle distribution function. 

Potential energy ([1]) can be written as 

UNiqi,...,qN) = ^ (t){.qjAk), (5) 

^<3<k<N 

where 

,/ A ( i\ , ui{q)+ui{q') 

Introduce a function h{q, q') by the relation 

exp{-/30(g,g')} = ^{l + %,g)}, (7) 



where 

1 



J exp{-p(f){q,q')}dqdq. (8) 



The canonical Gibbs distribution (12]) takes the form 



DNiqi,...,qN) = Qn n + (9) 



with 



Qn = Yi + h{(lj^(lk)]dqi- ■ -dqN- (10) 

l<j<fc<iV 

From ([9]) and ffTOl) it follows that statistical properties of system under consideration 
are completely determined by the specifying single function of two configurational vari- 
ables h{q, q') and two external parameters N and V . Gibbs distribution ([9]) has the same 
form as for system without external field. Therefore we may suppose that results of article 
[1] can be proved for system considered in this article too. 

In section [2] we state a mathematical problem for our molecular system and formulate 
conditions for existence and uniqueness of its solution. In section [3] feasibility of these 
conditions for considered physical system are proved. In section H] an expression for 
function h{q,q') in terms of F2{q,q') is calculated. In section [5] an expression for the 
canonical distribution in terms of reduced reduced two-particle distribution function is 
produced. In section E] it is shown that a set of first reduced distribution function and pair 
correlation function contains the whole information as well as two-particle distribution 
function and canonical Gibbs distribution. 

2 Mathematical formulation of problem 

Let us introduce function f{q,q') by the relation 

F2{q,q') = ^^^^^,[l + f{q,q% (H) 

Both the function f{q,q') and the function h{q,q') satisfy the conditions 

J f{q, q') dqdq' = 0, j h{q, q') dqdq' = 0. (12) 

From expressions Q, and (ITT!) it follows that 

1 + /(gi,92) = 7— / 11 [1 + h{qj,qk)]dq3- ■ -dq^. (13) 
^^■^ l<J<t<N 

This relation defines the transformation {h f} and can be considered as a nonlinear 
equation relative to h{q,q'). If there exists a solution h{q,q'] [/]) of this equation then 
a function Dj^^qi, . . . ,gjv) becomes an operator function of /. It means that both the 
canonical Gibbs distribution D^^ and all reduced distribution functions Fi are expressed 
in terms of a single reduced distribution function F2. Thus we have to prove that equation 



f|T3|) has a unique solution and therefore the transformation {h —>■ /} has inverse one 
Multiplying equation f[T^ by QnV ^ and using ([TU]) we rewrite it in the form 

[1 + / (91,52)] yTjv / n [^ + h{Pj,Pk)]dpi---dpN 

i<j<k<N 

Yl [1 + h{qj, qk)]dq3 ■ ■ ■ rfg^ = 0. (14) 

^<j<k<N 

The left-hand side of (IT^ is a polynomial operator of degree A/" = A^(A^ — l)/2 relative 
to h and of degree one relative to /. Denote this operator by J-'{h, /). Equation (fT^ can 
be written in a symbolic form 

^(/i,/) = 0. (15) 
To solve this equation it is necessary to specify an additional condition 

jr(/,(o) j(o)) ^ 0, (16) 

where f^^\q,q') and h^^\q,q') are assigned functions. 

We can easily determine these functions for our physical system. If the external field 
and interaction between particles are absent, i.e., potentials ui and U2 are constant, then 
the function h{q, q') vanishes. Under this condition Qn = Q^^ = , Di^{qi, . . . , g^r) = 
Df = V-^, and /(g, q') = 0. Therefore we can take = and /(°) = in (fT6l). 

Equation (|T5l) and additional condition ( JT6l) form a problem on implicit function. In 
functional analysis there is a number of theorems on implicit function for operators of 
various smoothness classes. We use the theorem for analytic operator in Banach space in 
the form given in the book [5]. 

Theorem. (On implicit function). Let J-'{h, f) he an analytic operator in D^{h!~^\ Ei) x 

Dp{f^^\ E) with values in E2. Let an operator B =^ —dJ-'{h^^\ f^^^)/dh have a bounded in- 
verse one. Then there are positive numbers pi andri such that a unique solution h = xif) 
of the equation J-'{h, /) = with the additional condition J-'{h^^\ f^^^) = exists in a solid 
sphere Dj.^{h^^\Ei). This solution is defined in a solid sphere Dp^{f^^\E), is analytic 
there, and satisfies the condition h^^^ = xif^^^)- 

Here ^^(^^o, ^) denotes a solid sphere of radius r in a neighborhood of the element Xq in a 
normalized space S, the symbol x denotes the Cartesian product of sets, dJ-'/dh is Frechet 
derivative |5l |6] of the operator J-', h^^^ and f^^^ are assigned elements of the respective 
spaces El and E. If the functions h{q, q') and /(g, q') satisfy this theorem conditions, the 
former is a single valued operator function of the later. 

To prove an existence and uniqueness of the solution of problem (1151) . flTB]) we have to 
show that conditions of the above theorem are satisfied. 

3 Proof of feasibility of the theorem conditions 

First we define spaces Ei, E, and E2 mentioned in the theorem for functions describing 
the physical system under consideration. 
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3.1 Functional spaces of problem 



Potentials Ui{q) and U2{q,q') are real symmetric functions. Suppose they are bounded 
below for almost all {q, q'} G V. All physically significant potentials possess this property. 
Under this condition integral ([8]) exists and h{q, q') is a real symmetric function bounded 
for almost all {g, g'} G V . 

A set of functions bounded nearly everywhere forms a complete linear normalized 
space (Banach space) with respect to the norm [HI E] 

= vraisup |/i(g,g')h (17) 

(<?,g')eV(2) 

where "vraisup" denotes an essential upper bound of the function on the indicated set 
and V*-^^ V yiV is Cartesian product of V by itself. It is called the space of essentially 
bounded functions and is denoted by Loo(V^^-'). In addition h{q,q') satisfies condition 
f|T2|) . The set of such functions is a subspace of Loo(V*^^^). It is easy to show that this 
subspace is a complete space relative to norm ( fTTl) . Therefore we can take the Banach 
space of symmetric essentially bounded functions satisfying condition ( fT2l) as Ei. 

Expression ( |T3l) for /(g, q') includes multiple integrals of different power combinations 
of h{q,q'). Any power of essentially bounded function are integrable with respect to 
arbitrary set of variables {qji, ■ ■ ■ ,qjr} over V [8\. Therefore all integrals in (fT^ are 
essentially bounded functions too. Arguing as above, we can show that the space E of 
functions f{q,q') coincides with Ei. Continuing in the same way we can show from fll4p 
and (1151) that E2 is the same space. Thus we define the spaces of the above theorem as 
E = Ei = E2 = Loo(V(2)) with property (fT2D . 

From ([7]) and ffTTj) for h{q, q') and f{q, q') it follows that / > — 1 and h > —1. Therefore 
we can take a manifold {/ > —1, h > —1} as a definition domain of the operator J-'{h, /). 
Since the left-hand side of f lT^ is a polynomial, the operator J-'{h, f) is analytical in 
this domain. As stated above the additional condition flTBl) is valid for h^^^ = f^^^ = 0. 
Thus any solid spheres of £"1, E with centers at h^^^ = 0, f^^^ = and radii r < 1, 
p < 1 respectively can be used as domains Dr{h^^\ Ei) and Dp{f^^\E) indicated in the 
theorem. 

Finally it is necessary to prove that the operator 

dcf dT{h, f) 
^ = dh- 

has a bounded inverse one. 

3.2 Properties of the operator B 

To find the inverse operator it is necessary to solve the equation Bh = y, where 
h & El, y & E2. The expression for Bh is a linear relative to h part in the left-hand side 
of relation f|T^ as / = 0. Let us introduce next notations for arbitrary function ^(g, g') 

l{q) = lpj e(g, q')dq', ? = ^ / ^(^' (19) 

Expanding products in f|T4|) and keeping linear summands we obtain 

{Bh){q, g')) = h{q, q') + {N - 2)[h{q) + h{q')]. (20) 



f = 



It's easy to estimate a norm of the operator B. Using definition f|T7|) we obtain 



\Bh\\ < {2N-3)\\h\\. 



From here we get an estimation 



1511 < 2A^-3. 



Thus the operator B is bounded. 

Using (1201) we can write a nonuniform equation Bh = / in the next form 

h{q,q') + {N - 2)[h{q) + h{q')] = f{q,q'). 

Taking into account the conditions (fT2l) we get a solution of this equation 

N-2 



hiq,q) = iB-'f){q,q) = fiq,q')- 



N 



[f{q) + f{q')]- 



(21) 
(22) 

(23) 
(24) 



From here it's easy to estimate a norm of the inverse operator B ^. Evaluating the norm 
of right-hand-side of (12^ we get 



IB-'II < 



\\B-'f\\ 





N- 


2 




< 












1 





3(A^-1) -2 



N-1 



< 3. 



(25) 



Therefore the operator B^^ exists and it is bounded. 

So all conditions of the above theorem are valid for our physical system. Hence there 
exists a unique solution h = xif) of problem ( ITSl) . ( |T6|) as a function of /. This solution 
defines an inverse transformation from the function / to the function h{q, q') = h{q, q'; [/]). 



4 Derivation of the inverse transformation h{f) 

To obtain the transformation f h we have to solve equation (|T^ relative to h{q,q'). 
Here we present less unwieldy derivation than in the work P]. At first define an auxiliary 
operator function g{h) by means of a relation 

Yl [l + h{qj,qk)]dq3---dqN = l + g{h). (26) 

This operator function is a polynomial of degree Af relative to h and depends on two 
configurational variables qi,q2- It can be written in the form 

g{h) = Y,9i{h), (27) 

1=1 

where gi{h) is a uniform operator of order / relative to h. Let us derive the expression for 
gi{h) from definition fl26l) . 

For the sake of abbreviation of subsequent calculations we introduce next notations. 
We will denote by number K & {!,••• , A/"} every ordered collection {j, k} from the set 
{!,..., A^}. Such one-to-one correspondence can be always made. A collection {g^, g^} 
is an element of manifold V*^^-*. We will denote this element by Xk- By definition put 
Xi = (gi,g2)- 




Expanding the product in fl26|) we obtain an expression 

9iih) = [dq3--- dq^ ^ KXk,) ■ ■ ■ KXk,) (28) 

for every / = !,... Introduced operators gi{h) as well as g{h) are symmetrical func- 
tions of two configurational variables: gi{h) = gi{q,q'; [h]) and g{h) = g{q,q'; [h]) . In 
contrast to h and / both g{q,q'; [h]) and gi{q,q'; [h]) don't satisfy condition (fT2|) except 
for gi{q,q'; [h]). First term of series (l27|l is gi{h) = Bh and satisfies to condition (fT2l) . 
Configurational integral (ITOl) takes the form 

A/" 

Q7v = V^^(l + ^) = ^^[l + $^^.]. (29) 

k=2 

Here we used notations ( |T9l) . 

Substituting definitions fl2Bl) and fl2Ul) into f[T^ we write it in the form 

(l + /)[l + ^(/i)]-[l + (7(/i)] = 0, (30) 

where / and g{h) are functions of Xi = (gi,g2)- But value g doesn't depend on configu- 
rational variables, it is a functional relative h. Substituting expansions fl271) and fl29l) here 
we reduce this equation to the form 

/[I + Y^m)] + Y.'^m - gi{h)] -Bh = 0. (31) 

1=2 1=2 

Here we have taken into account that gi{h) = Bh and 'gi{h) = 0. 
For subsequent calculation we need multilinear operators 

Gi{yi,...,yi) = I dqs-'-dqN ^ yi{XK,) ■ ■ ■ yiiXxi). (32) 

These operators are linear with respect to any functional argument We can consider 
the operator functions gi{h) as generated by these multilinear operators Gi 

gi{h) = Gi{h,...,h). (33) 

Operators Gi{yi, ■ ■ ■ ,yi) are functions of configurational variables {gi, 52} = Xi. In gen- 
eral these functions aren't symmetrical relative to {qi,q2)- But this isn't important since 
under substituting of these operator functions into equation fl3T|) symmetric property will 
be hold automatically. In the result we can rewrite equation fl3Tl) as 

M _ Af _ 

h = B-'f[l + J^GiiK ...,h)] + B-' . . . , /i) - Giih, ...,h)]. (34) 

1=2 1=2 

We will search a solution of this equation in the form of power series 

h = Y,h,{f), (35) 

k=l 
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where h^-lf) are uniform operators of order k relative to /. At the same time they 
are functions of configurational variables Xj. Substituting (1551) into (1511) and taken into 
account linearity of Gi{yi, . . . ,yi) with respect to any argument yi we obtain 

oo Af CO oo 

k=l 1=2 ji=l 31=1 

M oo oo 

+B~' E E ■ • ■ YP^^^n^ • • • ' ^n) - Gi{h,,,. • • , /.J]. (36) 

1=2 h=l ji=l 

Transform sums over ji , . . . , ji as follows 



oo oo 



ji=l j;=l k=l jiH hi!=fc 

Then relation (136|) takes the form 

oo Af OD 

hif) = B-'f + E E E ^^(^^- • • • ' ''n)! 

k=l 1=2 k=l jiH yii=k 

J\f oo 

+B-'Y.Y1 E [Gi{h,,,...,h,,)-Gi{h,,,...,h,,)]. (37) 

/=2 fc=i j^ + -+ji=k 

Double sum Y^^=2 ^T=i transformed as follows 

Af oo Af oo N-l k oo M 

1=2 k=l 1=2 k=l 1=2 k=N k=2 1=2 k=M 1=2 

or as follows 

Af oo Af Af Af oo Af k oo Af 

EE = EE+E E =EE+ E E- 

1=2 k=l 1=2 k=l 1=2 k=Ar+l k=2 1=2 k=Af+l 1=2 

Substituting these expressions into (1571) we get the relation 

oo Af k—1 

hkif) = B-'f + E E E ^'(^^- • • • ' 

k=l k=3 1=2 hi;=fc-l 

oo Af 

+B-' J2 H E G,ihn^...,h,,)f 

k=X+l 1=2 ji+---+ji=k-l 
Af k _ 

+B-'J2Y1 E [Gi{h,,,...,h,,)-G^{h,,,...,h,,)] 

k=2 1=2 ji+-+ji=k 
oo Af 

+B-' E E E mhn,---,h,,)-G^{h,„...,h,^)]. (38) 

k=Ar+l 1=2 ji + -+ji=k 
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Here in the first two sums of tlie riglit-liand side we cliange summation variable k to k + 1. 

In tliis relation all sums with respect to k contain expressions of order k relative to /. 
Putting terms of the same order being equal in accordance with the theorem on uniqueness 
of analytical operators [9J we obtain the next recurrent system for the functions hk{f) 



hif) = B-'f, (39) 
h2{f) = B-'[G2{h, h) - G2{h, h^)], (40) 

k-l k 

hkif) = B-'Y. E ^'(^^1' • • • ' ^n)f + 

1=2 ji+---+ji=k-l 1=2 

■ J2 [^iihn.---,h,,)-Giih,,,...,h,,)l 3<k<J^, (41) 

iiH l-jl=k 

_ AT 

h,{f) = B-'Y1 E • • • ' ^^•^)/ + E 

1 = 2 jl + ...+j;=fc-l 1 = 2 

■ J2 Pi{h,,,---,h,,)-Gi{h,,,...,h,,)l k>M+l. (42) 

iiH ^3i=k 

All terms of series ( l35l) are calculated from this system. So the solution of the equation 
(ITSll is founded. It satisfies additional condition (|T6l) . Convergence of the series (|35ll 
with hk{f) being the solutions of system ( |39l) - (H2l) is proved by Cauchy-Goursat method 
presented in the book 

As soon as hk are expressed in terms of / we can get the canonical distribution ([9]) in 
terms of F2 since / and F2 are uniquely bounded by relation (|TT|) . 



5 Calculation procedure for the canonical distribu- 
tion in terms of / 

Since canonical distribution ([9]) is a ratio of two polynomials with respect to we see that 
Z^AT is an analytical operator function of h. We have just proved that h is an analytical 
operator function of /. Therefore Z^at is an analytical operator function of / and it can 
be expanded into an absolutely convergent series relative to / 



Dm = V-'' 



k=l 



where fk{f) is a uniform operator of order k transforming function f{q,q') to function 
(fkiQi, ■ ■ ■ ,(lN][f])- Taking into account definitions of the reduced distribution functions 
dl]) and the function /(g, q') f[TT|) we can get relations for ipk 

jdqs--- dqN^iiqi, . . . , gjv; [/]) = f{qi, ^2), (44) 



yN- 



Jdq3---dqN^k{qi,---,qN][f]) =0, k = 2,3,.... (45) 
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Below we construct a procedure for calculation of functions (pi{qi, ■ ■ ■ , Qn] [/]) in terms of 
/• 

We introduce a nonlinear operator function X{h) by the relation 

J] [l + h{q,,qk)] = l + X{h). (46) 

This operator function is a polynomial of degree Af relative to h. It can be written in the 
form 

M 

X{h) = Y,Mh), (47) 



1=1 



where Xk{h) are defined by relations 



AKgi,...,giv;[/^])= Yl KXk,) ■ ■ -hiXK,). (48) 

l<Kx<-<Ki<N 

Introduce also multilinear operators 

AKyi,...,yz) = yi{XK,)---m{XK,). (49) 

l<Ki<---<Ki<J\[ 

It's evident that 

Xi{h)=ki{h,...,h). (50) 

The operators introduced here are connected with the operators giji)^ gk{h), and 
Gk{hi, . . . , hk) by the relations 

dq^--- dqNX{qi, . . . , g^; [/i]) = g{qi, q^] [h]), (51) 



yN-2 
1 



yN- 



dqs--- rfgjvAfc(gi, . . . , g^v; [/i]) = gkiqu 92; [h]), (52) 
1 



In particular for k = 1 



dqs+i ■ ■ ■ dq^Akiqi, . . . ,qN;[hi, . . . , hk]) 

= Gk{qi,...,qs;[hi,...,hk]). (53) 



dq3---dqNAi{qi,...,qN; [h]) 



yN-2 

Giiqu q2] [h]) = giiqi, q2] [h]) = {Bh){qi, ga)- (54) 



Taking into account the expression fl29|) for Qn we can write 



. V-'^lt^. (55) 

l + g{h) 



Comparing it with we get the relation 

HI 

+m 



1 + O /l 



fc=l 

in 



where h is the operator function of / calculated in previous section. Using here the 
expressions for A(/i), 'g{h) and h{f) we can transform the right-hand side of ( l56l) to series 
with respect to / and thus obtain expressions for v?fc(/). But less awkward transformations 
are obtained if we construct a recurrent system for ifk{f)- 

Multiplying (!56l) by 1 + g{h) and using (133|) and (1501) we obtain 

AT _ oo M ^ _ 

{1 + Y^GkiK h)} J2 viU) = E ^'^(^' ^) - (5^) 

fc=2 Z=l k=l k=2 

Substitution of the expansion (l35l) here gives 

A/" oo oo oo 

k=2 ji=l jk=l 1=1 
A/' oo oo AT oo oo 

= E E ■ ■ ■ E Mhn,-, h,,) - E E ■ ■ ■ E ^^■^)- (58) 

fc=i ii=i ifc=i fc=2ji=i jk=i 

Further calculation is carried out in the same way as in the previous section. We won't 
make it and write a recurrent system for V2fc(/) straight away 



<^i(/)_= Ai(/ii), (59) 
^2(/) = ^l{h2) + ^2{hl. - h,), (60) 

k _ 

Mf) = Mhk)+Y, E {Hhn,-,hj,)-Gi{hj„...,hj,)} 

1=2 ji+-+ji=k 
k _ 

-E E Gz_i(/i,„...,/i,,_J^,,(/), 3<A:<Ar, (61) 

'=3 iiH hji=k 

N _ 
ifkif) = Ai{hk) + Y1 E {Mhn,...,hj^)-Gi{h,,,...,h,,)} 

1=2 jj+-+ji=k 
A/'+l _ 

-E E GUhn,-.h^rWM, k>Af+l. (62) 

'=3 iiH \-ji=k 

In this relations we have to use the expressions for hr{f) derived from the recurrent 
system (I39D-(!12]). 

For example an expression for (fi{f) is 

N 

^i{qi,---,qN;[f])= Yl fi^vqk)-{N-2)Y,J{q,). (63) 

l<j<k<N i=l 

This expression coincides nominally with (y9i(gi, . . . , g^v; [/]) derived in the paper [4J. How- 
ever in this formula the function /(g, q') depends on both the interaction potential U2{q, q') 
and the external field ui{q) whereas in [4] the latter is absent. It is easy to show that 
expressions ( !59|) - (l62!) and (163|1 satisfy conditions ^M . (H5|l . 

1 1 



6 Set of irreducible functions containing the whole 
information about system under consideration 



The program presented here can't be reahzed for the reduced one-particle distribution 
function because not all the theorem conditions are held in this case. Namely the operator 
Bi = dFi/dh{q,q') has a nontrivial space of zeroes. This space consists of all functions 
h{q, q') G Ei satisfying a condition 

j h{q,q')dq' = 0. (64) 

Such operator B might not have an inverse one. And the reduced one-particle distribution 
function doesn't contain the whole information about system under consideration. 

At the same time there is a set of irreducible functions which is equivalent to the 
reduced two-particle distribution function. This set includes the reduced one-particle 
distribution function and a pair correlation function 



— mww — ■ 

The function Fi (g) can't be expressed in terms of q') and vice versa. 

Taken together they form the set which describes the system under consideration com- 
pletely because there is a one-to-one correspondence -^2(5, q')< — >{Fi{q), x(g, q')}. Really 
on the one hand F2{q,q') is expressed in terms of mentioned set of functions by relation 
fl^Sjl . On the other hand Fi{q) is expressed as an integral of F2{q,q') with respect to q'. 
And x(g, q')} is defined by relation (165!) in which Fi{q) is presented in terms of F2{q, q'). 

So the function f{q,q') can be expressed in terms of the set {Fi{q), g')}. There- 
fore the canonical Gibbs distribution can be expressed in terms of this set of irreducible 
functions with the help of the relations (|^ - (|^ in which functions hi{q,q'; [/]) must 
be expressed in terms of {Fi{q) and ^{q,q')}- Therefore this set of irreducible functions 
contains the whole information about the system under consideration. 



7 Conclusion 

Using the theorem on implicit functions in this article it is shown that the reduced dis- 
tribution function of order two plays a specific role for the canonical ensemble of N 
particles with two-body interaction and external field. The canonical Gibbs distribution 
D^iqi ■ ■ ■ , (?Ar) can be expressed uniquely in terms of this function ^2(5, q'). From here 
we easily conclude that there is a one-to-one correspondence between these two functions. 
This means that the reduced distribution function F2 contains information about the 
system under consideration as much as the whole canonical distribution D^. Reduced 
distribution functions of all orders can be expressed in terms of this single function F2. 

The one-particle reduced distribution functions don't satisfy the theorem conditions. 
So it is impossible to express the canonical distribution in terms of this function Fi. To 
all appearance it contains not all information about the system under consideration. 

The set of irreducible functions {Fi(g), x(g, q')} as well as -^2(5, q') contains the whole 
information about the system under consideration. 
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Considered theorem provides sufficient conditions for existence and uniqueness of in- 
verse transformation {/ —>■ h}. Results obtained here are vahd in some neighbourhood 
of /t(°) = 0, f^^^ = 0. The question about size of this neighbourhood demands special 
investigation. 
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